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Abstract 

We give a review of the state of the art with regard to the theory of scale 
functions for spectrally negative Levy processes. From this we introduce a 
general method for generating new families of scale functions. Using this 
method we introduce a new family of scale functions belonging to the Gaussian 
Tempered Stable Convolution (GTSC) class. We give particular emphasis to 
special cases as well as cross-referencing their analytical behaviour against 
known general considerations. 
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1. Spectrally negative Levy processes and scale functions 

Let X = {X t : t > 0} be a Levy process defined on a filtered probability space (Q,!F, F,P), 
where \Tt ■ t > 0} is the filtration generated by X satisfying the usual conditions. For x E R 
denote by P x the law of X when it is started at x and write simply Po = P. Accordingly we 
shall write K x and E for the associated expectation operators. In this paper we shall assume 
throughout that X is spectrally negative meaning here that it has no positive jumps and that 
it is not the negative of a subordinator. It is well known that the latter allows us to talk 
about the Laplace exponent ip(0) := \ogE,[e eXl ] for $1(6) > where in particular we have the 
Levy-Khintchine representation 

i)(6) = -a6+ ]-(T 2 e 2 + [ (e 9x - 1 - x ei {x> - 1} )U(dx) (1) 

where a € M, a > is the Gaussian coefficient and II is a measure concentrated on (— oo,0) 
satisfying j,^ j (1 hx 2 )W(dx) < oo. The, so-called, Levy triple (a, a, II) completely characterises 
the process X. 

For later reference we also introduce the function $ : [0, oo) — > [0, oo) as the right inverse of 
ip on (0, oo) so that for all q > 

$(<?) = sup{(9 > : ip(6) = q}. (2) 
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Note that it is straightforward to show that ip is a strictly convex function which is zero at the 
origin and tends to infinity at infinity and hence there are at most two solutions of the equation 

m = q. 

Suppose now we define the stopping times for each i£R 

r+ = inf{t > : X t > x} and r" = inf{f > : X t < x}. (3) 

A fluctuation identity with a long history concerns the probability that X exits an interval [0, a] 
(where a > 0) into (o, oo) before exiting into (— oo, 0) when issued at x g [0, a]. In particular it 
is known that 

where x g (— oo, a], q > and the function : M — > [0, oo) is defined up to a multiplicative 

constant as follows. On (— oo, 0) we have W^\x) = and otherwise is a continuous 

function (right continuous at the origin) with Laplace transform 

i 

e - ex W {q) {x)dx = —— for ^(9) > (5) 

tp(0) - q 



The functions {W^ : q > 0} arc known as scale functions and for convenience and consistency 
with existing literature we write W in place of W^°\ Identity (4) exemplifies the relation between 
scale functions for q = and the classical ruin problem. Indeed setting q — we have that 
^x(tq < T a ) — 1 — W(x)/W{a) and so assuming that ip'(0+) > and taking limits as a | oo it 
is possible to deduce that W(oo) _1 = -0'(O+) and hence 

P x (r " < oo) - 1 - V/(0+)VK(x). (6) 

It is in this context of ruin theory that scale functions make their earliest appearance in the works 
of Zolotarev (1964), Takacs (1966) and then later either explicitly or implicitly in the work of 
Emery (1973) and Korolyuk (1974, 1975), Suprun (1976) and Rogers (1990). The real value of 
scale functions as a class with which one may express a whole range of fluctuation identities for 
spectrally negative Levy processes became apparent in the work of Chaumont (1994, 1996) of 
Bertoin (1996, 1997) and an ensemble of subsequent articles; see for example Lambert (2000), 
Avram ct al. (2004), Pistorius (2003, 2004, 2005, 2007), Kyprianou and Palmowski (2006), Doncy 
and Kyprianou (2006) and Doney (1991, 2005a, 2007). Moreover with the advent of these new 
fluctuation identities and a better understanding of the analytical properties of the function 
came the possibility of revisiting and solving a number of classical and modern problems from 
applied probability, but now with the underlying source of randomness being a general spectrally 
negative Levy processes. For example, in the theory of optimal stopping Avram et al. (2002, 
2004); Alili and Kyprianou (2005) and Kyprianou (2006), in the theory of optimal control Avram 
et al. (2007) and Loeffen (2007), in the theory of queuing and storage models Dube et al. (2004) 
and Bekker et al. (2007), in the theory of branching processes Bingham (1976) and Lambert 
(2007), in the theory of insurance risk and ruin Chiu and Yin (2005), Kluppelberg et al. (2004), 
Kluppelberg and Kyprianou (2006), Doney and Kyprianou (2006), Kyprianou and Palmowski 
(2007), Rcnaud and Zhou (2007), Loeffen (2007) and Kyprianou et al. (2008), in the theory of 
credit risk Hilberink and Rogers (2002) and Kyprianou and Surya (2007) and in the theory of 
fragmentation Krell (2007). 
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Although scale functions are now firmly embedded within the theory of spectrally negative 
Levy processes and their applications, and although there is a reasonable understanding of how 
they behave analytically (see for example the summary in Kyprianou (2006)), one of their main 
failings from a practical point of view is that there are limited number of examples which are 
available 'off-the-shelf for modeling purposes. 

The original purpose of this paper was to introduce a previously unknown family of scale 
functions. However, in doing so, we came across a number of additional examples which were 
seemingly less well known in the literature or are themselves new and concurrent with this article. 
We shall therefore initially spend some time documenting known examples of scale functions. It 
seemed then sensible to also document known general analytical properties of scale functions and 
this we have also done in Section 3. 

Following that we shall expose a new family of explicit examples of scale functions as well 
specifying in detail the associated Levy processes. By 'explicit' we mean that, as in the below, 
known examples, we are able to give exact analytical expressions for the scale functions. By 
'specifying in detail' we mean that we shall describe the Levy triple (a, a, II) of the associated 
spectrally negative Levy process and in particular it will turn out that a density for the Levy 
measure can be obtained and we are able to identify it as a known function. Moreover, it turns 
out that the associated Levy processes discussed may also be described as the independent sums 
of other, well known Levy processes. 

2. Known examples of scale functions 

At the time of writing, and to the best of the authors' knowledge, there are essentially six main 
examples of explicit classes; three of which deal with the case of a compound Poisson process 
with negative jumps and positive drift. 

1. Brownian motion with drift. The associated Laplace exponent is written ip(6) = a 2 9 2 / '2 + fi8 , 

where a > and \i £ R then, for x > 0, 

W {q) (x) = 2 e -A"^ 2 sin h(^ ^2qa 2 + Ji) (7) 
V 2q<7 2 + /i cr 

for q > 0, where in the case that q = fi = the above expression is to be taken in the 
limiting sense. 

2. Spectrally negative stable processes with stability parameter (3 £ (1,2). When ?/>(#) = 9 13 

with (3 £ (1, 2), we have for x > 

W^ q \x) =/3x - 1 E' f3!l (qx ) (8) 

for q > 0. where Ep t i(z) = J2k>o z k /T(l + (3k) is the Mittag-Lcfflcr function. 

There also exists an expression for the scale function of the aforementioned /3-stable process 
but now with a strictly positive drift (but only for the case q = 0) implicitly in the paper 
of Furrer (1998). Indeed, by considering the expression there for the ruin probability one 
may extract the following, 

W(x) = - c (l-E^ 1 (-cx^ 1 )) (9) 

for a: > 0. Note that by taking c | one obtains an expression which agrees with (8), 
namely x ~ l /T{(3). 
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3. Generally speaking it is possible to deal with a strictly positive linear drift minus a compound 

Poisson process process with positive jumps with a rational Laplace transform. The 
calculations are however rather extensive. See for example Mordecki and Lewis (2005) 
where the scale function can be extracted from their expression for the ruin probability. In 
the special case that the drift is denoted c and jumps arrive at rate A and are exponentially 
distributed with parameter /i such that c— Xj [i > 0, that is to say ip(9) = c9 — A(l — //(// + 
O)^ 1 ), we have for the case q = and x > 0, 

W(x) = - fl + — ^— (1 - e<"- c " lA > x )y (10) 
c y C[i — X J 

The expression for the scale function in the case that q > with exponentially distributed 
jumps, and indeed for the slightly more general case where an independent Gaussian 
component is added in, can in principle be extracted from Kou and Wang (2003). 

4. This example is taken from Asmussen (2000) and again comes from the case of a known 

ruin probability giving a scale function for the case q = and positive drift. Consider a 
spectrally negative compound Poisson process whose jumps are exactly of size a G (0, oo), 
with arrival rate A > and with positive drift c > such that c — Xa > 0. In that case, 
V>(0) =c6- A(l - er a9 ) and for x > we have 

W{x) = _ £ e-^-^ c - ( - J (an-xr (11) 

n=l ' ^ ' 

where [x/a\ is the integer part of x/a. 

5. The next example comes from Abate and Whitt (1999) who themselves refer to older work 

of Boxma and Cohen (1998). 

Consider a unit-rate linear drift minus a compound Poisson process of rate A > with 
jumps whose distribution F has Laplace transform given by 

/•oo n 

/ e~ 6x F{dx) = 1 = (12) 

for 3?(6!) > The latter Laplace transform corresponds to a random variable whose mean 
is equal to > and whose tail distribution takes the form 

F(x,oo) = (2x + l)rj{x) - 2y^, x > (13) 
when \i = 1 and otherwise when [i ^ 1 

F(x, oo) = (jzr^) ~ ^( x ^ 2 ))' x ^ ( 14 ) 

where rj{x) — e x : erfc( v / a;) . It is assumed that the underlying Levy process drifts to oo. 
Since E(Xi) = 1 — X/fi the latter assumption is tantamount to A//i < 1. 

Said another way, we are interested in a spectrally negative Levy process with Laplace 
exponent given by 

<P{6) = 8 p ^. (15) 

(p + Ve)(i + Vo) 
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For the scale function with q = it is known that for x > 

w{x) = T=hfr i 1 W^V^ 11 ^ - ^(^ 2 i))) ■ (16) 



and 



6. The final two examples, both scale functions only for the case q = 0, appeared very recently 
in the theory of positive self-similar Markov processes, see Chaumont et al. (2007). See 
also Chaumont and Caballero (2006) for the origin of the underling Levy processes. The 
first example is the scale function which belongs to a spectrally negative Levy process with 
no Gaussian component, whose Levy measure takes the form 



e (0-i)y 

n{dy)= (ev~ ip+i dy ' y< ° (18) 



where /3 € (1,2) and whose Laplace exponent takes the form 

for $1(0) > 0. Note that tp'(0+) < and hence the process drifts to — oo. In that case it 
was found that for x > 

W(x) = (1 - e~ x ) - 1 e x . (20) 

The second example is the scale function associated with the aforementioned Levy process 
when conditioned to drift to oo. It follows that there is still no Gaussian component and 
the Levy measure takes the form 

U{dy) = (e«- dy ' V < ° (21) 
and the associated Laplace exponent is given by 

m = mm 

for $t(0) > 0. The scale function is then given for x > by 

W(x) = (1 - e~ x f-\ (22) 



3. Known analytic properties of scale functions. 

Although at the time of writing further examples other than those above are lacking, there are 
a collection of general properties known for scale functions, mostly concerning their behaviour at 
and oo. For later reference in this text and to give credibility to some we review them briefly 
here. As usual, (a, a, LI) denotes the Levy triple of a general spectrally negative Levy process. 
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Smoothness The following facts are taken from Lambert (2000), Chan and Kyprianou (2007), 
Kyprianou et al. (2008) and Doney (2005b). It is known that if X has paths of bounded variation 
then, for all q > 0, W^|(o,c») € C 1 (0, oo) if and only if II has no atoms. In the case that X has 
paths of unbounded variation, it is known that, for all q > 0, W^|(o,oo) e C 1 (0, oo). Moreover 
if a > then C 1 (0,oo) may be replaced by C 2 (0,oo). It was also noted by Renming Song (see 
the remarks in Chan and Kyprianou (2007)) that if II(— oo, — x) is completely monotone then 

^ (9) |(0,oo) G C°°(0,Oo). 

Concavity and convexity It was shown in Loeffen (2007) that the latter assumption that IT(— oo, —x) 
is completely monotone also implies that W^'(x) is convex for q > 0. Note in particular, the 
latter implies that there exists an a * > such that is concave on (0, a*) and convex on 

(a*,oo). In the case that ?//(0+) > and q = the argument in Chan and Kyprianou (2007) 
shows that a* — oo and W is necessarily concave. In Kyprianou et al. (2008) it is shown that if 
II(— oo, — x) has a density on (0, oo) which is non-increasing and log-convex then for each q > 0, 
the scale function (x) and its first derivative are convex beyond some finite value of x. 

Continuity at the origin For all q > 0, 

, . I if a > or f, , n .(— x)H(dx) = oo 

W (q) (0+) = {_!.. n , J( 7^ 0)V . ' ^L.t , ^ (23) 
v ' I c 1 it a = and J^_ 1 Q ^(— x)U(dx) < oo, v ' 

where c = —a — J^_ 1 ^ xX\{dx). 

Derivative at the origin For all q > 0, 

2/cr 2 if a > 

W iq) '{{)+) = { oo if cr = and n(-oo, 0) = oo (24) 

( (? + n(-oo,0))/c 2 if a = and n(-oo,0) < oo. 



Behaviour at oo /or g = As x \ oo we have 

nrt ^ / 1 /^'(0+) if ^'(0+) > 

~ | e *(0)x/^/ ($(0)) if v y (0+) < o. (25) 

When E(-Xi) = 0a number of different asymptotic behaviours may occur. For example, if writing 
tp(6) = 9(f)(9) (which is possible by the Wiener-Hopf factorization), we have </>'(0+) < oo then 
W{x) - x/0'(O+) as x t oo. 

Behaviour at oo for q > As x f oo we have 

ff(«»(i)~e* ( '»7f'(*((i)) (26) 
and thus there is asymptotic exponential growth. 

4. Methodology for new examples 

For any given spectrally negative Levy process, the scale functions are intimately connected to 
descending ladder height process and this forms the key to constructing new examples. Therefore 
we shall briefly review the connection between scale functions and the descending ladder height 
process before describing the common methodology that leads to the new examples of scale 
functions. For a more detailed account of this connection, the reader is referred to the books of 
Bertoin (1996), Kyprianou (2006) or Doney (2007). 
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It is straightforward to show that the process X — X_ := {X t — X_ t : t > 0}, where X_ t := 
inf s < t X s , is a strong Markov process with state space [0,oo). Following standard theory of 
Markov local times (cf. Chapter IV of Bertoin (1996)), it is possible to construct a local time 
at zero for X — X_ which we henceforth refer to as L — {L t : t > 0}. Its inverse process, 
L- 1 := {L^ 1 : t > 0} where L t _1 = inf{s > : L s > t}, is a (possibly killed) subordinator. 
Sampling X at L^ 1 we recover the points of minima of X. If we define H t — X L -i when 

L^ 1 < oo with H t = oo otherwise, then it is known that the process H = {H t : t > 0} is 
a (possibly killed) subordinator. The latter is known as the descending ladder height process. 
Moreover, if T is the Levy measure of H then 

/■OO 

T(x, oo) = e^°> l e-*( ^II(-oo, -u)du for x > 0, (27) 

J X 

see for example Vigon (2002a). Further, the subordinator has a drift component if and only if 
g > in which case the drift is necessarily equal to er 2 /2. The killing rate of H is given by the 
constant E(Xi) V 0. 

The starting point for the relationship between the descending ladder height process and scale 
functions is given by the Wiener-Hopf factorization. In 'Laplace form' for spectrally negative 
Levy processes this takes the appearance 

#) = («-$(0M«) (28) 

where 

4>{6) = -logE(e- eHl ), il>{0) =logE(e eXl ), (29) 

and $1(6) > 0. With this form of the Wicncr-Hopf factorization in mind, we appeal principally 
to two techniques. 

1. We choose the process H, or equivalently <p(9), so that the Laplace inversion of (5) may 
be performed. In some cases, when the process does not drift to — oo (or equivalently 
f/>'(0+) > 0), it can be worked to ones advantage that, after an integration by parts, one 
also has that 

r°° 1 

I e ~ 6xw ^ = W) m 

for Sft(0) > 0. 



2. We choose the process H to be such that its semigroup ¥(H t G dx) is known in explicit form 
and then make use of the following identity (see for example Bertoin (1996) or Kyprianou 
(2006)), 

f-OC 

I dt ■ P(H t e dx) = W(dx), x>0, (31) 
Jo 

whenever X does not drift to — oo. 



Naturally, forcing a choice of descending Ladder height process, or equivalently requires one 
to know that a spectrally negative Levy process, X, exists whose Laplace exponent respects 
the factorization (28). The next Theorem provides the necessary justification. Indeed, to some 
extent it shows how to construct a Levy process with a given descending ladder height process 
as well as a prescribed ascending ladder height process. For spectrally negative Levy processes, 
the ascending ladder height process is a (killed) linear unit drift and has only a single parameter, 
namely $(0), the killing rate. Hence the construction in the below theorem offers a parameter (p 
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which plays the role of <&(0). Since both ascending and descending ladder height process cannot 
both be killed one also sees the parameter restriction 

(pK = (32) 

where k is the killing rate of the descending ladder height process. 

Theorem 1. Suppose that H is a subordinator, killed at rate n>0, with Levy measure which is 
absolutely continuous with non-increasing density and drift (. Suppose further that ip > is given 
such that ipK = 0. Then there exists a spectrally negative Levy process X , henceforth referred 
to as the 'parent process', such that for all x > 0, P(t+ < oo) = e~ vx and whose descending 
ladder height process is precisely the process H. The Levy triple (a, a, II) of the parent process is 
uniguely identified as follows. The Gaussian coefficient is given by a = \/2C The Levy measure 
is given by 

II(-oo,-a:) = <^T(x,oo) + ^^. (33) 

dx 

Finally 

a= xU(dx) - k (34) 

J (-00,-1) 

if ip = and otherwise when ip > we can establish the value of a from the eguation 

a<p=\a 2 <p 2 + [ (e** - 1 - x<pl {x> _ 1} )IL(dx). (35) 

Z J(-oo,0) 

In all cases, the Laplace exponent of the parent process is also given by 

m = (8- <p)4>{6) (36) 

for6>0 where (f)(9) = - logE( e - 0Hl ). 

Proof. The proof is reasonably self evident given the preceding account of the ladder height 
process. One needs only the additional information that if X is any spectrally negative Levy 
process with Levy measure II then 



E(Xi) = -a+ xll(dx). (37) 

J (-00,-1) 

Moreover when $(0) > then necessarily the descending ladder height process has no killing 
and when the descending ladder height process is killed then $(0) = 0. Further, for all x > 
P(t+ <oo) ^ e -*(°K □ 



The idea of working 'backwards' through the Wiener-Hopf factorization as we have done above 
can also be found in Bertoin et al. (2004) and Vigon (2002b). Note that it is more practical to 
describe the parent process in terms of the triple (a, II, ip) than the triple (a, a, II) and we shall 
frequently do this in the sequel. It is also worth making an observation for later reference 
concerning the path variation of the process X for a given a descending ladder height process H. 

Corollary 1. Given a killed subordinator H satisfying the conditions of the previous Theorem, 
(i) the parent process has paths of unbounded variation if and only if T(0, oo) = oo or ( > 0, 
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(ii) if Y(0, oo) = A < oo then the parent process necessarily decomposes in the form 

X t = (K + X- (p)t + y^CB t - S t (38) 

where B = {B t : t > 0} is a Brownian motion, S = {S t : t > 0} is an independent driftless 
subordinator with Levy measure v satisfying v(x,oo) = II(— oo, — x). 

Proof. The path variation of X follows directly from (33) and the fact that a = Also 
using (33), the Laplace exponent of the decomposition (38) can be computed as follows with the 
help of an integration by parts; 

p oc p oc rl'Y* 

(k + A — (p)0 + CO 2 -<f>0 e' 8x T(x, oo)dx -6 e- 9x — {x)dx 

Jo Jo dx 

C°° dT C°° dT 

= (k + T(0,oo) - (<p)6 + (O 2 -if J {I - e- 0x )^-{x)dx - 9 J e - 0x — {x)dx 

- (0-tp) ^ + 6( + J™(l-e- Bx )^-(x)dx^ . 

This agrees with the Laplace exponent ip{6) = (6 — (fi)(f>(9) of the parent process constructed in 
Theorem 1. □ 



5. The Gaussian tempered stable convolution class 

In this paper we introduce a new family of spectrally negative Levy processes from which our 
new examples of scale functions will emerge. We call them Gaussian tempered stable convolution, 
GTSC for short. When there is no Gaussian part we call the distribution a tempered stable 
convolution and write TSC for short. 

5.1. The tempered stable ladder process 

The starting point for the construction is a tempered stable subordinator plus a linear drift, 
possibly killed, that will play the role of the descending ladder height process for the GTSC parent 
process. Some references for tempered stable distributions and tempered stable Levy processes 
are Steutel and van Harn (2004), Cont and Tankov (2003) and Schoutens (2003). Note, that 
several names and origins, and many different parameterisations are used for the tempered stable 
distributions and processes. 

In our parametrization, the tempered stable subordinator involves three parameters, the 
stability parameter a < 1, the tempering parameter 7 > 0, and the scaling parameter c > 0. 
When a < it is required that 7 > 0. Furthermore we might add a linear drift with rate ( > 
to the process, and possibly kill the process at rate k > 0. 

The Laplace exponent of the tempered stable subordinator, and henceforth the descending 
ladder height process is thus taken to be 1 

4>{6) = n + Ce + cT{-a)( 1 a -( 1 + e) a ), R(0)>- 7 , (39) 

and the associated Levy measure is given by 

T{dx) = cx-^e'^dx (x > 0). (40) 

For < a < 1 the process has infinite activity. For a = the expression (39) is to be understood 
in a limiting sense, i.e., <j>(0) = n + (8 — clog (7/(7 + 9)), and the process is simply a (killed) 



1 Note that we use T(z) as a meromorphic function with simple poles at the non-positive integers. By the 
functional equation of the gamma function we have r(-a) = -a _1 r(l - a), and T(-a) < for a 6 (0, 1). 
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gamma subordinator (with drift). If a < the process is a compound Poisson process with 
intensity parameter cT(— aYf a and gamma distributed jumps with —a degrees of freedom and 
exponential parameter 7 > 0. 

5.2. The associated parent process 

We may now invoke Theorem 1 to construct the associated GTSC process. This introduces 
another parameter ip > 0. To be able to apply the theorem, we need a decreasing Levy density, 
and thus we have to restrict the scaling parameter to — 1 < a < 1. The theorem tells us that the 
Laplace exponent of the parent process is 

1>(0) = (k - yQ6 + (9 2 + c(6 - <p)T(-a)( 7 a - (7 + 6) a ). (41) 

for 31(9) > 7, and obviously for a = we understand the above expression in the limiting sense so 
that ip(6) = (K + (ip)8 + (9 2 — c(6 + ip) log (7/7 + 9). It is important to recall here and throughout 
the remainder of the paper that rnp = 0. The corresponding Levy measure given by 

U(dx) = c ^^K ^dx + c-^-zU^cto (42) 

for x < 0. This indicates that the jump part is the result of the independent sum of two spectrally 
negative tempered stable processes with stability parameters a and a + 1. We also note from 
Theorem 1 that a = \/2Cj indicating the presence of a Gaussian component. This motivates the 
choice of terminology Gaussian Tempered Stable Convolution. 

If < a < 1 the jump component is the sum of an infinite activity negative tempered 
stable subordinator and an independent spectrally negative tempered stable process with infinite 
variation. If a = the jump component is the sum of a spectrally negative infinite variation 
tempered stable process with stability parameter 1 and exponential parameter 7 and the negative 
of a gamma subordinator with exponential parameter 7. If — 1 < a < the jump part of the 
parent process is the independent sum tempered stable subordinator with stability parameter 
1+a and exponential parameter 7, and an independent negative compound Poisson subordinator 
with jumps from a gamma distribution with —a degrees of freedom and exponential parameter 
7. In the extreme case a = — 1, the parent process has negative jumps which are compound 
Poisson and exponentially distributed with parameter 7. 

6. Evaluating GTSC scale and q-scale functions 
6.1. The case with rational nonzero stability parameter 

Theorem 2. Suppose a = m/n with m e Z, n e Z, < \m\ < n, and gcd(|m|, n) = 1. Let us 
consider the GTSC process with Laplace exponent 

V>(0) = (9 - tp) [k + (9 + cT(-a) (7" - (7 + 8) a )] . (43) 

Let m + = max(m, 0) and m_ = max(— m, 0). Then the polynomial 

f q (z) = [z n - 7 - tp) [(k + (7 + cT(-a) 1 a )z m - + ( z n + m - - c r(-a)z m +)] - qz m - (44) 

has at least one real root. Let I denote the number of distinct roots, n , . . . , rg the roots, arranged 
such that r\ is the largest real root, and fii, . . . , fj,£ their multiplicities. Let Akj for k — 
j = 0, . . . , jJLk — 1 denote the coefficients in the partial fraction decomposition 
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Then we have 

*(«) = r" - 7 (46) 

and 

W^) (a .) = e -7- £ £ A fc ^x^+ 1 )/"- 1 M 3 \(r fe2; ^). (47) 
fe=l j=0 ■ ? ' 

Proof. We have 

/,(*) = M*" - 7) - 9) (48) 
for z e Z with | arg(z)| < ir/n, and conversely, 



1 (0 + 7)" 



W-9 /,((* + 7)») 



(49) 



for G C\ (— oo, —7]) such that the denominator on the left hand side does not vanish. We know 
from the theory of scale functions that 9 = <&(g) is a real root of tp(6) — q = 0, and there is no 
larger real root of that equation. Thus z = (&(q) +7)" is a real root of f g (z) — and there is 
no larger real root of that equation. This shows (46). 

The left hand side of (49) is the Laplace transform of W^ q \x), and we know it is an analytic 
function for $t(0) > Using the partial fraction decomposition we get 



/ e -¥«'(# = yy , (50) 

Jo ~k 



and this equation is valid for any q > and 3?(6*) > Q(q). We recognize the Laplace transform of 
exponentially tilted derivatives of Mittag-Leffler functions on the right hand side, see for example 
(Jacob, 2005, Prop.7.1.9, p.359), and obtain (47). □ 



Generally speaking this is rather an awkward formula to work with analytically. The main 
strength of the expression lies with it being a simple matter to program into a package such as 
MATLAB or Mathematica. It is instructive to revisit the (known) analytical properties of W(x) 
and W( q \x) listed in section 3 and check them directly for our explicit expression as this will 
yield in most cases more detail. 

Smoothness: Since the Mittag-Leffler functions are entire functions, we see that W^(x), which 
is a priori defined for x € (0, 00), admits an analytic continuation to x € C \ (—00, 0], and thus 
is C°° on (0, 00). 

Behaviour at zero: We obtain the behaviour at zero from the power series expansion of the 
derivatives of the Mittag-Leffler function. Exploiting algebraic relations of the and Akj several 
terms cancel and we obtain for £ > 

W<«>'(aO~i (i-»0). (51) 
A similar argument applies for £ = 0. If < a < 1 we get 

T a T a—1 
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In the case that — 1 < a < and £ = 0. we get again by cancellation of terms and limiting 
behaviour at zero of Mittag-Lcffler functions 

W {q \x) l — -, wW{x)~-. CX " 1 (i-»0). (53) 

V ; K + cj a T(-a) V ; (k + C7°T(-q;)) 2 ' V J K ' 

Note, another way to establish these last two results is by Karamata's Taubcrian Theorem and 
the Monotone Density Theorem respectively from the asymptotics of ip(6) as 9 — > oo. This 
approach shows, that the asymptotics above also hold for irrational a e (— 1, 1). 

Behaviour at infinity: The behaviour at infinity can be obtained from the asymptotics of the 
Mittag-Lefflcr functions at infinity. The dominating contribution comes from the term with n. 
It is useful to recall that f g (z) = ip(z n — j) — q and thus f' q {z) = nz n ~ 1 ip'(z n — 7). Let us consider 
first q = 0. Note, that we have ip'(0+) = k - tp(( + C7 a_1 r(l - a)). Suppose ip'(0+) > 0, which 
implies that ip = 0. Then n = 7™ is a single root of fo(z). After some elementary simplifications 
we obtain 

W(x) ~ - (x -» 00) (54) 

in agreement with the theory. Suppose now ip'(0+) < 0. Then n = (</? + 7)" is a simple root 
of fo(z). As above, using the asymptotics of the Mittag-Leffler functions at infinity, we see that 
the term with r\ dominates all other terms. After elementary simplifications we obtain 

£ <fiX 

W(x) — (x -» 00). (55) 

V ; K + C<^ + cr(-a)( 7 « - ( 7 + (^)«) y J y ' 

Consider the third case, ^'(0+) = 0. This can only happen when k — and ip = 0. Then n = 7™ 
is a double root of fo{z). The asymptotics show, that the contribution from n dominates and 
we obtain 

^)- C + cr 4(l-a) (^°°)- (56) 
Finally in the case g > 0, little more can be said than (26). 

6.2. The case with general nonzero stability parameter a E (— 1, 1) 

When we take g = and ( = and a e (—1, 1)\{0} (without the restriction of begin a rational 
number) one may obtain much cleaner expressions for the scale function than the formulation in 
Theorem 2. 

Theorem 3. Suppose a e (—1, 1) \ {0} and consider a TSC process (without Gaussian compo- 
nent) and Laplace exponent 

m = {6-<p)[ K + cT(-a) ( 7 Q - (7 + 0) a )} ■ (57) 

Then If < a < 1 i/ien 

j/ — 1 < a < 0, i/ien 
W(z) 



K + cr(-a)7 a 

*(-<*- r.-.^^-.^f j » ,59) 



e 



(K + cT(-a)j a ) 2 J a "\K + cT(-a)r 
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Proof. This follows from the known Laplace transforms 

f°° 1 
J e- ex x a - 1 E a , a (\x a )dx=-^-j (60) 

and 

J™ e-^X^x-^E.^X^x-^dx = - 1, (61) 

valid for a > resp. a < together with the the well-known rules for Laplace transforms 
concerning, primitives and tilting. □ 



With these closed-form expressions, the explicit asymptotics discussed after Theorem 2 are 
immediate by inspection. Note in the particular case that k = — cL(— a)j a the expression for 
the scale function in the case < a < 1 reduces simply to 

W ^ = ~ I" ^TTV^e-^dy. (62) 
K J 1 (a) 

Note that although the Laplace inversions involved in the proof of Theorem 3 are straightfor- 
ward, it is not clear that the resulting expressions for W are scale functions without the presence 
of Theorem 1. 

6.3. The case with a — i; inverse Gaussian descending ladder height process 

We will assume now that £ = 0, tp = 0, k = 0. Substituting a i— » 1/2, c i— » 5/V2n, 7 1— > 7 2 /2 
reveals that the tempered stable subordinator we have been working with is in fact the familiar 
inverse Gaussian IG(S, 7)-subordinator when a = 1/2. It turns out for this case that we also get 
some cleaner expressions than the one given in Theorem 2. 

When the descending ladder height process is the IG(5,j) subordinator, the jump part of 
the parent process is a superposition of a negative IG subordinator and a spectrally negative 
tempered stable process with stability parameter 3/2. 

Following the general approach in Theorem 2 we relate ip(6) to a polynomial, which, in the 
IG parametrization of this section, is f q (z) = ip(z 2 — j 2 /2) — q. The first simplification is, that 
f q (z) is a polynomial of degree 4 if £ > 0, and degree 3, if ( = 0. Thus completely explicit, 
elementary (and somewhat lengthy) expressions of its roots can be given in terms of radicals 
by the formulae of Ferrari resp. Cardano. The second simplification concerns the Mittag-Lefflcr 
functions. As a = 1/2 the Mittag-Lefflcr function simplifies to an expression involving the more 
familiar (complementary) error function. 

We find that fo{z) has a positive double a negative single root, and after some elementary 
simplifications, 

1 



W(x) = — [(1 + l 2 x) crfc (^-7 V /a72j + 7VW^ e " 57 x - l\ ■ (63) 

Actually, in this special case the result can be verified faster by calculating the Laplace transform 
of (63). As a side remark, we note, that W"{x) = — .T~ 3 / 2 e~27 x /(2^/2tt8) and thus we see, that 
W is indeed concave, as we already know from the theory. Let us now consider q > and put 



q = i<5 7 3 . 



(64) 
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If < q < qo then we have three simple real roots ri,r2,rs, that can be obtained by Cardano's 
formula 2 , and 



W {q \x) = e~?' 



r 3 e r 3 x 



W) Crfc (-^) + ^) -fc(-r 2 ^) + erfc(-^) 



(65) 



If q > qo we have a simple real root r\ and two complex conjugate simple roots and r^. The 
same formulas hold as in the previous case. Alternatively, we could in the complex case write 



W (q \x) = e -^ 2x 



— ^- erfc(-nv^) + 23? < — - — - erfc(-r 2 \/i) 



/i(r 2 ) 



(66) 



If q = qo then we have a simple positive and a double negative root, and after some easy 
simplifications, 



W (qo \x) 



1 



36<S 7 



67 y| e -^ + l 5e l^erfc(-|y| 
-e-^ 2;r (15 + 27 2 a;) erfc 



7 p 
3V2 



(67) 



6.4. The stable case 

Because of their definition via their Laplace transform, it is immediate that scale functions 
are continuous in the parameters of the parent process. Taking 7 J. in the Laplace exponent 
of the tempered stable ladder height process shows us that the parent process is an a + 1 stable 
process when k = <p = C = and c = — 1/T(— a). Taking limits as 7 J. in the expression for the 
scale function in the first part of Theorem 3 gives us the known result W(x) — x a /T(a + 1). 

6.5. The case a — 0: gamma descending ladder height process 

When a = we can only treat the special case with q — 0, k = 0, ( = and <p = 0. 

Theorem 4. 

/>oo 

W(x) = / P{ct,-yx)dt, (68) 
Jo 

where P(a, x) is the regularized lower incomplete gamma function 

Proof. We may appeal to the second method in (31) using resolvents to note that the scale 
function density satisfies 

W'{x) = -x^e^^-log^x)) (69) 



2 Explicitly 



= _j_ 2v^5 7 2 3/- 8 ^3 7 3 +27 ^2 g+3 y3 v ^ 7 ^4 g 2_ 16g 5^ 

3v/5 3 ^/-8«3 7 3 + 27 1 52 9+ 3%/3v / 27<S4< 7 2_ 164 5 ; ^ 3V2S 

V2(l+is/3)5j 2 (l-iV3) Y/-8i 3 7 3 +27« 2 q+3v / 3V27i 4 9 2 -165 5 7 3 <7 

3v/ 2 3 ^/-8« 3 7 3 +27 1 52 9+ 3v^V /27 ' s4 9 2 -16« 5: ^ 

_ 2 v / 2(l-iv / 3)*7 2 (l+iV3) Y / -8i 3 7 3 +27i 2 q+3\/3v'27« 4 q 2 -16<S 5 7 3 q 

3v/5 3 ^/-8« 3 7 3 +275 2 < 7 +3v^ V / 2 7 ' s4 9 2 -16<S s 7^9 6 ^ 
If q = <57 3 /2 those formulae simplify, as one root is then zero. 



15 



where 

roo -dx 

<PV) = / j^rdx, (70) 
Jo T ( x ) 

which is the Laplace transform of the reciprocal gamma function. 3 As X has infinite variation, 
we have W(0+) = 0. Integrating (69) yields the result. □ 

6.6. Remarks on the cases a — —1/2 and a — — 1 

The case a = — 1/2 is again covered in Theorem 2 above, but several simplifications occur. 
The jump part of the ladder process is in this case a compound Poisson process with jumps from 
the xi-distribution. Again, we can avoid Mittag-Leffler functions and use the complimentary 
error function instead. With regard to the roots of the polynomial equation wc may say the 
following. If C > the polynomial f q (z) has degree 5 and its roots typically cannot be expressed 
in terms of radicals. If ( = the polynomial is of degree 3 and Cardano's formula can be used. 
If in addition n = and ip = then we have r\ = y^y a double and r 2 = —y/i a single root. 
We can proceed as is Section 6.3 to obtain some simple closed form expressions for W(x) and 

When a — —1 the jump part of the parent process is a compound Poisson process with negative 
exponential jumps, and the ladder process is a compound Poisson process with exponential jumps. 
The Laplace transform of the scale functions are rational, and the scale functions can be obtained 
directly from the partial fraction decomposition. The results are implicitly contained in Mordecki 
and Lewis (2005) and Kou and Wang (2003) and we exclude the calculations here. 

7. Numerical illustrations 

Figure 1 contains graphs of the GTSC scale function W(x) for stability parameters a = 1/4, 
1/3, 1/2, 2/3, 3/4 in six cases. Figure 2 contains graphs of the GTSC (/-scale functions W^ q \x) 
with q = 1, but otherwise for the same parameters. Below we introduce six cases considered in 
terms of the classification of the parent process. 

• Case A, k = 0, <p = 0, £ = 0, c = 1, 7 = 1: The parent process is oscillating, has no 
diffusion part and infinite variation jumps. The ladder process is an infinite activity pure 
jump subordinator, has no linear drift, and is not killed. 

• Case B, k = 1, ip = 0, £ = 0, c = 1, 7 = 1: The parent process drifts to +00, has no 
diffusion part and infinite variation jumps. The ladder process is an infinite activity pure 
jump subordinator killed at rate k = 1. 

• Case C, k = 0, <p = 1, £ = 0, c = 1, 7 = 1: The parent process drifts to —00, has no 
diffusion part and infinite variation jumps. The ladder process is an infinite activity pure 
jump subordinator. 

• Case D, k = 0, <p = 0, £ = 1, c = 1, 7 = 1: The parent process is oscillating, has no linear 
drift, but a diffusion part and infinite variation jumps. The ladder process is an infinite 
activity pure jump subordinator plus a linear drift. 

• Case E, k = 1, ip = 0, £ = 1, c = 1, 7 = 1: The parent process drifts to +00, has 
a Gaussian component and infinite variation jumps. The ladder process is an infinite 
activity subordinator plus linear drift, killed at unit rate. 



3 The value at = is called the Franscn-Robertson constant and denoted by <p(0) = F. A lot of analytical 
and numerical material on tp{6) can be found in Franscn and Wriggc (1984). That paper contains also several 
references to related work going back to Paley, Wiener, Hardy and Ramanujan. Note, that the reciprocal gamma 
distribution in Fransen and Wriggc (1984) is not to be confused with the inverse gamma distribution, which is 
often called reciprocal gamma distribution as well. 
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• Case F, k = 0, <p = 1, £ = 1, c = 1, 7 = 1: The parent process drifts to —00, has a 
diffusion part and infinite variation jumps. The ladder process is an infinite activity pure 
jump subordinator plus a linear drift. 

Now let us discuss the graphs in view of the theoretical properties listed in Section 3. 
Smoothness: All graphs look smooth, as the theory predicts. 

Concavity and convexity: For q = we know in cases A, B, D and E the graph is predicted 
to be concave. In cases C and F the graph is predicted to be convex-concave. Indeed this is the 
appearance in Figure 1. 

For q = 1 we know, that all graphs are concave-convex. We observe this behaviour in Figure 2. 
In some cases, it is necessary to inspect the shape of the graph more closely in the neighbourhood 
of the origin in order to see concavity. 

Behaviour at zero: For all cases we observe W(0) = in Figure 1 and W^ q \0) = in Figure 2, 
in agreement with formula (23). 

For the cases A-C (no diffusion part) we observe W'(0) = +00, for the cases D-E (nonzero 
diffusion part) we observe W'(0) = 1 in Figure 1. We observe the same behaviour, IfW'(O) = DO 
resp. W^'(0) = 1, in Figure 2. This is in agreement with formula (24). 

Behaviour at infinity: Firstly let us consider the case q = 0. In cases A and D we have 
y/(0+) = and thus, according to (56), asymptotic linear growth as x — > 00. In cases B and E 
we have n > and hence -0'(O+) > so that W(x) -^1/Kasnoo. In cases C and F we have 
ip > and hence ip'(0) < and so there is exponential growth of W(x) as x — > 00 according to 
the second case in formula (25). This corresponds to what is observed in the graphs in Figure 1. 
(The behaviour for x — > 00 becomes more prominent when plotting < x < 20.) 

For the case q > 0, all graphs in Figure 2 appear exhibit exponential growth as x — > 00, in 
agreement with formula (26). 

8. Additional Remarks 

8.1. Extensions and ramifications 

Our methods apply also to the integrated q-scale functions, see (Kyprianou, 2006, Theorem 8.1, 
p.214) 

Z^{x) = 1 + q [ W {q) {y)dy, (71) 
Jo 

for the Gaussian tempered stable convolution class with rational stability parameter. 
Also we could treat more complicated Gaussian tempered stable convolutions with 

n 

4>(6) = k + (9 + J2 Cin-aiW - ( 7i + B) ai ) (72) 

»=1 

with a\, . . . , a n rational numbers in [—1,1). 

The key observation is simply, that any Laplace transform, that is of the form f(9 a ) with 
/ a rational function and a some number, can be inverted in terms of the partial fraction 
decomposition of / using Mittag-Leffler functions and its derivatives. 

In other very recent work, Kyprianou and Rivero (2007) have looked at families of scale 
functions which are the result of choosing the descending ladder height process to have a Laplace 
exponent which is a special or complete Bernstein function. 
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Figure 2: Scale functions W {q) (x) for the GTSC class with q = 1 and < a < 1 
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8.2. Numerical inversion of the Laplace transform 

While our formulas provide closed form expressions for rational a in terms of partial fraction 
decompositions and Mittag-Lefner functions, those expression can be rather lengthy, involve the 
zeros of polynomials of high degree, and, experience shows them to be numerically problematic 
when a — m/n unless n is quite small. 

Alternatively one may always consider numerical inversion of the Laplace transform of W^ q \ 
which is extremely simple for the GTSC class. Moreover this approach works also for irrational 
a. We use the Bromwich inversion integral for the Laplace transform, 

r+zoo 

W M {x) = J_ f J^L (73) 
K ' 2m J t/>(0)-q V ' 

r — ioo 

where r > &(q). Let us consider the integrand for \0\ — > oo with $1(0) > r. The term e 0x is 

bounded. If C > we have ip(6) ~ C<? 2 , if C = and < a < 1 we have tp(0) cT(-a)9 a+1 

and (73) is a proper Lebesgue integral. Otherwise it must be interpreted as a principal value 
integral. Some remarks how to handle this situation and how to improve the numerics and a 
piece of Mathematica code, that evaluates the scale function by (straightforward, and numerically 
unsophisticated) integration, are in the extended report. For further material on the numerical 
evaluation of scale function see also (Surya, 2007, Chapter7, p. 11 Iff). 

8.3. Examples of applications 

Let us briefly conclude by re-enforcing the genuine importance of establishing explicit examples 
of scale functions from the point of view of modelling by giving three classical examples from the 
theory of applied probability. It will suffice to consider the role of 0-scale functions. We shall do 
so in the light of some of the explicit examples above. 

Example 1. Recent literature suggests that a more modern approach to the theory of ruin 
should involve considering, in place of the classical Cramer-Lundberg process, a spectrally 
negative Levy process which drifts to oo for the risk process. As noted earlier, for any such 
given Levy process with Laplace exponent if), the quantity 1 — ip'(0+)W(x) is the probability of 
ruin when the Levy process is issued from x at time 0. 
Consider the following parent process 

X t = ( K + \)t- Sf) - S* t (2) ,t>0 (74) 

where {S^ : t > 0}, i = 1,2, are subordinators whose Levy measures, i/W, i = 1,2, are identified 
by 

v {1 \dx) = X^-^x^e-^dx and v^(dx) = A^^y(l - v)x v - 2 e^ x dx (75) 

for x > 0. 

Following the interpretation of Levy risk processes in Kliippelberg and Kyprianou (2006) and 
Huzak et al. (2004) we may think of the parent process as the result of premiums collected at 
rate X + k to offset the result of two competing claim processes. The first has claims arriving at a 
compound Poisson rate A which are gamma distributed and the second has an infinite intensity 
of small claims whose aggregate behaviour is similar to that of a stable subordinator with index 
1 — v. Moreover, these claims are interlaced with larger claims of finite intensity, the tail of whose 
distribution are similar to those of gamma distributed jumps. 
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Recalling the discussion in Section 5.2 we have from the second part of Theorem 3 that the 
probability of ruin from an initial reserve x > is conveniently given by 

l -T^—-T^I X y-'e-^E^ipYvndy (76) 

A + K A + K J 

Example 2. If instead we consider the parent process in the previous example reflected in its 
supremum, then we are looking at the workload of a so called M/IT/1 queue (cf Dube et al. 
(2004)). That is to say a queue in which work is processed with a constant rate, say fi, and 
which arrives in packets according to a Poisson point process such that a job of size x > arrives 
in the interval (t, t + dt) with probability H(dx)dt + o(dt) where / °°(1 A x)U(dx) < oo. For the 
case at hand \x = A + n and H(dx) = u^\dx) + u^ 2 \dx) and there is again the possibility of 
interpreting the incoming work as the result of two competing mechanisms. In this case, the 
stationary distribution of the workload is given by one minus the ruin probability given above. 

Example 3. In the very recent work of Loeffen (2007), the classical dividend control problem 
of De Finetti's dividend problem has been considered in the context of the driving risk process 
being a general spectrally negative Levy process denoted as usual X. Here the objective is to 
find the optimal strategy and value to the following control problem 

v(x)=sapE x ^ e' qt dL^ (77) 

where q > 0, a™ = inf{i > : X t - < 0} and L n = {L^ : t > 0} is the process of dividends 
paid out associated with the strategy 7r and the supremum is taken over all strategies tt such 
that 17 is a non-decreasing, left-continuous adapted process which starts at zero and for which 
ruin cannot occur by a dividend pay-out. 

In Loeffen (2007), the author proves the remarkable fact that, under the assumption that 
the dual of the underlying Levy process has a completely monotone density, for each q > 0, 
the associated g-scale function, W^ q \ has a first derivative which is strictly convex on (0, oo). 
Moreover, the point a* := inf{a > : W^'(a) < W^'(x) for all x > 0} is the level of the 
barrier which characterizes the optimal strategy. The latter being that dividends are paid out in 
such a way that the aggregate process has the dynamics of the underlying Levy process reflected 
at the barrier a* . The value function of this strategy is given by 



W {q) {x) 



for < x < a 



C - a + WW (a) iorx > a 

If we are to take any of the parent processes mentioned in this paper as the underlying spectrally 
negative Levy process, then it is immediately obvious from (42) that they all have the property 
that their dual has a Levy density which is completely monotone. In such a case one may 
observe the value a* graphically (see for example Figure 2) as well as compute it numerically to 
a reasonable degree of exactness with the help of software such as Mathematica. 
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